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Abstract 

We discuss a framework relying on both perturbative and non-perturbative lattice 
computations which will be able to reliably determine the parameters of the EW phase 
transition. A motivation for the use of 3d effective theory in the lattice simulations, 
rather than the complete 4d one, is provided. We introduce and compute on the 
2-loop level a number of gauge-invariant order parameters - condensates, which can 
be measured with high accuracy in MC simulations. The relation between MS and 
lattice condensates is found, together with the relation between lattice couplings and 
continuum parameters (the constant physics curves). These relations are exact in the 
continuum limit. 
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1 Introduction 



The infrared problem in the thermodynamics of Yang-Mills fields prevents a purely 
perturbative study of the electroweak phase transition. It is well known that the 
effective potential for the scalar field cannot be computed perturbatively for small 
(j), in the symmetric phase. Accordingly, all the quantities characterizing the phase 
transition region, Tc, metastability and tunnelling temperatures, (j)c, latent heat and 
interface tension are not computable in perturbation theory. At the same time, some 
properties of the high energy particle excitations in broken and unbroken phases allow 
perturbative treatment. Moreover, the broken phase of the electroweak theory is weakly 
coupled, provided the vacuum expectation value of the Higgs field is large enough. Thus 
it is essential to apply both perturbative and non-perturbative methods to determine 
the nature of the electroweak (EW) phase transition. 

These general remarks suggest that probably the best way to determine the param- 
eters of the electroweak phase transition relevant for cosmology with sufficient accu- 
racy should consist in two steps, separating the perturbative and the non-perturbative 
physics. At the first step one goes as far as possible with analytical perturbative calcu- 
lations, which can simplify greatly the underlying 4d EW theory. The second consists in 
numerical lattice Monte Carlo simulations of non-perturbative physics. This program 
was initiated in ref. [1] and consists in the following ingredients: 

(i) Dimensional reduction of 4d high temperature theory. This step provides an 
effective 3d theory containing one (or even two) essential mass scales less than the 
underlying 4d theory. The scales of the 4d theory are the temperature T, the Debye 
screening mass rriD ~ gT, and the infrared scale gl = g^T , while those of the reduced 
theory are mj) and g\ (or just g\ if Aq - the temporal component of the gauge field - is 
integrated out). This part of the analysis has been carried out in [2] on the 2-loop level; 
the 2-loop contribution is essential both for numerical and conceptual reasons. The 3d 
theory is an excellent approximation to the 4d world for realistic Higgs masses. At the 
same time, this theory is much easier to analyze with perturbative or non-perturbative 
methods. The main reason is its super-renormalizable character. In contrast to the 4d 
theory, where ultraviolet divergences exist in any order of perturbation theory, in 3d 
only 1- and 2-loop graphs are divergent. This makes the scaling behaviour, essential for 
relating lattice and continuum parameters, much simpler than that in 4d, and allows 
one to make computations with high accuracy. 

(ii) A relation between lattice couplings and couplings of the 3d theory, which is exact 
in the continuum limit. With this relation, and using (i), one can relate any of the 
lattice observables to physical ones for given Higgs and W masses. The existence of an 
exact relation is a consequence of the super-renormalizable character of the 3d gauge- 
Higgs system, where 3-loop or higher order terms do not contribute to renormalization. 

(iii) The study of different observables on the lattice in the broken phase (below the 
critical temperature) and a comparison of them with 2-loop perturbative predictions 
^ . This allows one to estimate higher (3-loop, etc.) order perturbative corrections to- 

^As we have shown in [2] the 2-loop level of perturbative computation is a minimal one for which 
the result is practically scale-independent. 
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gather with an estimate of finite a (lattice spacing) effects. This part of the framework 
is quite important, since it can provide a justification for the use of perturbative meth- 
ods in the broken phase for concrete values of the Higgs mass and the temperature. 
If perturbation theory works with sufficient accuracy, then in some sense at least one 
half of the problems associated with the electroweak phase transition - the properties 
of the broken phase at Tc - can be dealt with analytical methods. 

(iv) The study of different observables in the symmetric phase (above Tc). Here the 
most interesting characteristics are correlation lengths of operators with different quan- 
tum numbers, in particular those that would be associated with Higgs and W bosons 
in the broken phase. The results of this study cannot be reproduced by perturbative 
methods, since the symmetric phase is in the strongly coupled confining phase. This 
part of the study is important since most of the mechanisms for electroweak baryoge- 
nesis are linked in one way or another with the properties of the symmetric phase. 

(v) As a combination of (iii) and (iv), the study of the system in the transition region 
and the determination of the parameters of the phase transition (such as its order, the 
critical temperature, the jump of the order parameter, etc.). 

This paper sets up points (ii) and partially (iii) of the above proposal, and adds 
several remarks to (i). These steps are necessary for the detailed study of the EW phase 
transition itself [3]. Some qualitative results on the EW phase transition has appeared 
in [1] (mostly for the light Higgs boson, niH = 35 GeV) and in [4] for niH = 80 GeV. 

The paper is organized as follows. In section 2 we provide further motivation of step 
(i), i.e. the use of a 3d effective theory [2, 4], [5] rather than the full 4d one [6]-[10]. In 
section 3 we summarize the results from the 2-loop perturbative analysis [2]. In sections 
4 and 5 we introduce a set of gauge-invariant observables (condensates), which can be 
computed with a simple manipulation from the 2-loop effective potential and can be 
measured on a lattice with high precision. A discussion of the renormalization of the 
ground-state energy Evac in the MS scheme is given in section 4. As a new perturbative 
result, section 5 contains a 2-loop perturbative computation of Evac and {(f)^(f)) in the 
MS scheme. 

The lattice action is given in Section 6, as well as a detailed discussion of relating 
lattice numbers to physical quantities, i.e. defining the constant physics curve. It is 
shown that the relation of the lattice to continuum requires the computation of three 
numbers, ?/, fj and rj. Two of these are computed in 2-loop lattice perturbation theory^. 
The third one, ?/, is related to the pure gauge sector, in which lattice perturbation theory 
is very complicated. We determine it in Section 7 by measuring {(f)^(f)) on the lattice 
deep in the broken phase and by comparing this with the 2-loop perturbative result. 
This comparison involves the constant physics curve and permits one to calibrate rj. 
The constant physics curves we found do not have any higher loop contributions^. 
Simultaneously, we use Monte Carlo data to estimate higher order corrections to the 

®In analogy to the derivation of the relation between Aqcd in lattice and continuum regularization 
schemes [11, 12] in QCD, the calculation amounts to computing the constant term besides the leading 
logarithmic term. However, now the logarithmic term arises only at the 2-loop level, which accounts 
for the greater complexity of the problem. 

'^In [1] the corresponding relation was obtained at the 1-loop level. 
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effective potential and the magnitude of finite- a effects. In Section 8 we construct the 
constant physics curves for a theory in which the Aq field - the former time component 
of the gauge field - is integrated out. Section 9 is a conclusion. Several technical 
points are covered in Appendices. In particular, we derive the 2-loop form of the 
gauge-invariant effective potential introduced in [13]. 

In this paper we frequently use the results of ref. [2]. All references to specific 
formulae from that paper are indicated by a I followed by the equation number. 

2 3d versus 4d 

In the following subsection we discuss the basic properties of dimensional reduction. 
We also introduce a criterion for a dimensional reduction different from that previously 
discussed in the literature and discuss the role of fermions. In the second subsection we 
discuss the quantitative advantages of the 3d formulation in the study of the electroweak 
phase transition. 

2.1 Why the 3d description is possible 

The idea of dimensional reduction [14]-[19] comes from a statement that equilibrium 
finite temperature field theory is equivalent to a Euclidean zero temperature field the- 
ory defined on a finite imaginary time interval l/T supplied with periodic boundary 
conditions for bosons and antiperiodic for fermions. Therefore, 4d finite temperature 
field theory is exactly equivalent to a 3d field theory at T = 0, but with an infinite 
number of fields. The 3d masses of bosons are ms = 2'KnT and those of fermions 
mp = TrT[2n + 1), n = 0, 1, .... 

The 3d effective bosonic action is defined as 

exp( — ^eff) = J D [heavy modes] exp( — 5), (1) 

where the functional integral is taken over all modes with masses ~ T. The effective 
action is a complicated functional of 3d bosonic fields. The use of it is, of course, 
equivalent to the full 4d theory. 

Suppose now that we are dealing with infrared physics only, so that the typical 
energy scale (which we denote as Q and discuss in more detail later) of the problem 
is much smaller than the temperature. Then all 3d fermionic modes and all bosonic 
modes besides the one corresponding to n = have masses larger than ttT and are 
very heavy in comparison with our scale. So, the integration over the heavy modes can 
be done perturbatively with an expansion parameter Q I{ttT). The effective 3d action 
can be represented in the form 

5eff = cVT' + j d\L^^{T) + Y.% (2) 

where Les{T) is a renormalizable 3d effective bosonic Lagrangian with temperature- 
dependent constants. On are operators of dimensionality n, suppressed by powers of 
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temperature, c is a number related to the number of degrees of freedom of the theory 
and V is the volume of the system. For example, the 3d gauge coupling is gl = g^T. 
The last step of dimensional reduction amounts to neglecting the terms Note that 
the constant part of the action is irrelevant for a discussion of phase transitions. 

Formally, the operators On are suppressed by the powers of temperature. It seems, 
therefore, that their contribution is negligible in the limit T — )• oo. This is, however, 
wrong since the 3d coupling and mass also grow without bounds in this limit. Thus 
On oc T", and the extra contributions to the 3d action do not vanish when T — )• oo. 
Thus, the Landsman [17] criterion of exact dimensional reduction 



where V^d is the effective potential computed in 4d and Vad is the potential computed 
with a renormalizable part of 3d action (with temperature-dependent parameters), does 
not hold. Instead, one has 



where the mi{T) are the relevant mass scales (inverse screening and correlation lengths) 
of the system. The question of validity of dimensional reduction thus becomes a prag- 
matic one: when is the r.h.s. of (4) negligible? At very high T all masses are propor- 
tional to coupling constants times T and the r.h.s. of (4) is some power of 4d coupling 
constants. These extra terms are small only when the dimensionless coupling constants 
(run to the scale of about TT, see below) of the underlying 4d theory are small. So, 
the criterion for dimensional reduction to work is thus the same as that for T = 
perturbation theory: g^ <^ 1 and A <C 1, A being a scalar self-coupling constant. For 
the electroweak theory (minimal standard model) this implies that mn should not be 
very large, say, mH<'imw 

However, we are mainly interested not in very large T but in T close to a possible 
phase transition temperature. Then the system contains two phases and may contain 
other mass scales than those set by the product of T and coupling constants. We will 
see below that the dimensional reduction is valid for the electroweak theory near Tc, 
provided that the Higgs mass is sufficiently large, say, mjf>30 GeV. 

One may wonder if we achieved anything going from 4d to 3d: dimensional reduction 
requires small coupling constants and works only when 4d perturbation theory works. 
The answer is: Yes, we do. The region of applicability of perturbation theory is different 
at zero and at high temperatures. The criterion for finite T perturbation theory to 
work is 



in hot QCD perturbation theory breaks down for distances Ij Q > l/[g'^T), etc.), while 
the criterion for T = perturbation theory to work, which is the same as the criterion 
of dimensional reduction to be valid, is g"^ <C 1. In other words, at high temperatures, 
we can find ourselves in a situation in which 4d finite T perturbation theory does not 
work {gl/Q = g^TjQ ^ 1), but dimensional reduction is possible (^f^ <C 1). As we 



lim {V^d - TVsd)/T' ^ 0, 

T—too 



(3) 



{V,d - TVsd)/T' = 0{mUT)/T') 



(4) 
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shall see, this is a typical situation in electroweak theory. After the reduction is done 
analytically, we can do some analytic perturbative or, most importantly, numerical 
non-perturbative analysis of the 3d theory. This may be easier and more precise than 
the analysis of the 4d one. 

Now, we can assert a different formal criterion for the validity of dimensional reduc- 
tion: take the limit T — )• oo, with 3d parameters fixed. Then the difference between 
4d and 3d computations always vanishes. This is not a physical limit since it implies a 
transition between different physical theories at zero temperature [g"^ = gl/T — )• 0). It 
is, however, quite helpful from a practical point of view, since it provides a definition 
of the formal expansion parameter, Q /T with Q ^ gl, which must be small for the di- 
mensional reduction be be valid. The important point is that in this limit the resulting 
3d theory has a non-trivial dynamics, which can be completely non-perturbative. 

A final comment concerns QCD at finite temperatures: one should not expect di- 
mensional reduction to be valid in the vicinity of the QCD phase transition, where 
the 4d gauge coupling is not small. However, at large enough temperatures, it may be 
used. 

2.2 Why a 3d description is more efficient 

In the previous subsection we argued that 3d effective theories can be used for the 
non-perturbative study of high temperature phase transitions in weakly coupled field 
theories. In this subsection we will discuss the qualitative and quantitative advantages 
of the 3d formulation. 

We start from the qualitative features. An important aspect of dimensional reduction 
concerns the role of fermions. Since the fermion masses (2n + l)TrT and bosonic non- 
static mode masses 2'KnT both get large simultaneously, both will give a contribution 
of the same order. Thus either both can be integrated out in the vicinity of Tc, i.e. 
the 4d— )-3d reduction is valid, or both have to be included and the 3d effective theory 
is not accurate*. Moreover, a realistic electroweak theory contains chiral fermions, 
which are impossible at present to put on the lattice. So, they must be integrated out 
analytically. In other words, 4d bosonic theory, which one can study on the lattice, is 
in any case an effective theory and not better than a dimensionally reduced 3d one. 

The second advantage is that complete 4d perturbative computations of the finite 
T effective potential are known only in the 1-loop approximation. The existing 2- 
loop computations [20]-[22] deal with high-temperature asymptotics only, which is 
equivalent to a 3d computation [2]. So, one is forced to compare the results of 4d 
lattice simulations with effectively 3d perturbative computations. Why then not carry 
out 3d simulations? 

The last, but not the least, advantage is the great simplification of perturbative 
computations in the 3d effective theory, in comparison with resummed perturbation 

*The former is the case for the electroweak theory for an interesting range of mjj. The latter 
holds for QCD, for which the gauge coupling is large at Tc. In fact, it is well known that the QCD 
transition depends sensitively on the fermion content and that, for pure SU(3) gauge theory, the 3d 
approximation starts to work somewhat above Tc, for T>1.5Tc. 
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theory in 4d. This allowed one to clarify the nature of the large 2-loop logarithmic 
corrections in the MS scheme and to introduce a method of summation of those based 
on the renormalization group [2]. 

We now turn to more quantitative aspects of the use of 3d effective theory. 



2.2.1 The quality of dimensional reduction in continuum 

We demonstrate our point by a study of different approximations to the quantity 
related to the 1-loop free energy of free bosons at high temperatures, entering the 1- 
loop evaluation of the effective potential and the computation of the effective scalar 
mass at high temperatures: 

where A;^ = + {2TrnTy. 

The high temperature expansion of (6) in continuum theory in the MS scheme reads: 



mT m\ ait . C(3) 



where = ^i^Te~'^ IT . The procedure of dimensional reduction takes into account 
the first three terms here. Namely, the first term changes the zero-temperature mass 
to the effective high-temperature mass, the third term fixes the scale at which the 4d 
coupling constants should be normalized, and the second term is purely 3-dimensional 
(it comes from the n = term in the sum). The 0{rn^) term and the higher-order 
terms are omitted in the effective 3d theory. They represent the contribution of higher 
dimensional operators to the effective action. 

By comparison of the second and the fourth terms in the expansion we see that the 
3d description is accurate within, say, a 1% level, provided that 



C(3) 



m 



3 



327r3 r3 



< 0.01, (8) 



i.e. m/T < 2.0. This condition is satisfied for all particles in the broken phase in 
the vicinity of the electroweak phase transition, if the Higgs mass is not very small, 
mn > 30 GeV. In fact, m/T is smaller than 0.5 for any realistic Higgs mass niH > 60 
GeV. This gives an accuracy better than 0.1%. In other words, the 3d description 
is a very good approximation to the real 4d high- temperature case, at least from 
the point of view of one-loop computations. The excellent convergence of the high 
temperature expansion has been mentioned in many papers on this subject (see, for 
example [23]-[25]). 

Of course, the small coefficients in front of higher powers of (m/T), found in one- 
loop approximation, do not necessarily reappear in 2 or higher loops. The accuracy of 
the dimensional reduction on the 2-loop level certainly deserves a careful investigation, 
which can be done on a basis of 2-loop computations in refs. [20, 21, 2, 22]. 
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2.2.2 Finite lattice spacing effects 

Next we wish to argue that the 4d lattice simulations of the electroweak phase transition 
can provide the same accuracy in the description of the EW phase transition as 3d 
simulations, only at the cost of a considerable increase of lattice volume. 

We start from a simple heuristic estimate. Let us take for simplicity a 4d lattice with 
equal lattice spacing a in temporal and spatial directions. The typical momentum of the 
particle in the plasma is about ST, and to have an adequate description of that on the 
lattice we must have a <C l/(3r). In contrast, in 3d simulations the only requirement 
is that the lattice spacing should be much smaller than a typical correlation length, 
a <C with m being a particle mass. For example, for the theory with Higgs mass 

niH = 80 GeV the ratio of W mass (defined as gv[T)/2 and measured on the lattice) 
to the temperature at the transition point is of the order of m/T ~ 0.2, so that the 3d 
requirement is about 3/0.2 = 15 times weaker than that for 4d simulations. Hence, in 
order to have the same finite a and finite volume effects, one should use 15^ ~ 3 • 10^ 
larger lattices for 4d simulations. 

This rough estimate can, in fact, be confirmed by an analytic study of finite a effects 
in the quantity I[m,T). For simplicity we take a spatially infinite lattice. Then, in a 
dimensionally reduced theory we have 

/dWd(m, T) = ^ + T/ - ^ log(^). (9) 

The function 7, being a 3d tadpole graph, has the following expansion (see Appendix 
A): 

I = - (am) - ^{amY + 0{amf], (10) 

47ra 

with S = 3.17591 and ^ = 0.15281 ^. Apart from counterterms, (9) differs from the 
continuum expression by two terms. First, it does not contain the 0{rn^ jT"^) term, 
dropped by the dimensional reduction procedure. As we argued above, this term is 
negligibly small for the realistic case of electroweak theory. Second, it contains terms 
proportional to a, which vanish in the continuum limit. We denote the magnitude of 
the finite a correction by p. To have, say, p = 10% (3%), finite a corrections to the 3d 
contribution to /s, the lattice should obey the condition l/(am) > ^jp ~ 1.5(5), which 
is quite easy to realize in practice. 

The 4d lattice approximation to (6) on an asymmetric lattice, which has a finite 
number of steps Nt = {aT)~^ in the temporal direction, is 

^t-l »^/a J3J^ I 

n=o J--/- i^'') + {-y E sm (^) + (-)2 sm (^) 
The large Nt (small a) expansion of this quantity reads: 

1 /I 94fi95\ 

W(m,r,iV,) = ^(^^) + + HNt)] 



^Contrary to the 4d case, where finite o corrections scale like o^, scaling violations start from linear 
in o terms in 3d. 
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+ T{h--— + ^-—) 12 
47ra Att 



2 /■fQ\ ™4 



m ,,,, _ C o m 



C(3) 



[log(iVO - 0.057 + f2{Nt)] + T^7^[l + fs{Nt)] + OK), 



where the functions fi[Nt), f2{Nt) and f3{Nt) vanish in the continuum limit iV^ — )• oo 
as 0{1/N^) and can be computed numerically (see Table 1). 





2 


3 


4 


5 


6 


7 


8 


9 


10 


fi 


0.200 


0.108 


0.056 


0.032 


0.021 


0.015 


0.011 


0.009 


0.007 


/2 


0.336 


0.144 


0.076 


0.047 


0.031 


0.023 


0.018 


0.014 


0.009 



Table 1: The values of functions fi and /2 for different numbers of sites Nt in the 
temporal direction. 

Let us first take for simplicity a lattice minimal subtraction scheme in which only 
the terms singular in a are removed by counterterms. The formal advantage of 4d is 
that the scaling violation is proportional to ~ and not to a, as in 3d. So, the 
functions fi and /2 decrease quite rapidly with increase of Nt. However, the numerical 
smallness of the coefficient ^ (reflecting the qualitative discussion of the different scales 
above), together with the smallness of the ratio m/T near the phase transition, makes 
the 3d case more advantageous than the 4d one. Indeed, to have the same level p 
approximation to 4d physics, one should have at least 

/i(iVO < ^P^. (13) 

For example, for p = 10% we have Nt > 5, and for p = 3% Nt > 8. Of course, this 
is not too difficult to achieve with computers available at present. However, at the 
same time, this requirement puts a rather strong limitation on the size of the lattice 
in the spatial direction. Indeed, for, say, p = 10% (3%) and m/T ~ 0.3 we have for 
the correlation lengths l/(am) ~ 16(55). In other words, to achieve in 4d the same 
accuracy of description of the EW phase transition as a 3d dimensionally reduced 
theory can provide, one should use about 10^ larger 3d volumes. In total, counting the 
number of sites in the 4th direction, the use of 3d theory provides a 3 to 4 orders of 
magnitude gain in computer time^°. 

The way out of this argument is the use of the temperature-dependent subtraction 
scheme, removing the finite contributions to /4diatt(m, T, Nt), proportional to fi and /2. 
This requires a careful study of the finite a behaviour of observables in 4d simulations 
at zero and non-zero temperatures^^. It is, however, highly non-trivial (but in principle 
possible [9]), since 4d theory contains at least three different mass scales (temperature, 
Debye screening mass, and 3d coupling g^T). 

^°If the computer time is not a problem, then the 4d simulations with, say, Nt = 8, may be 
superior to 3d simulations, provided the 2- and higher-loop dimensional reduction corrections come 
with coefficients so large that these corrections are larger that 4d finite o effects. 

^^We are grateful to I. Montvay for a discussion on this point. 
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3 Survey of continuum results 



It this section we summarize the continuum formulae relating 3d and 4d theories, and 
present a number of 2-loop relations in 3d. The starting point is the 4d Lagrangian of 
SU(2) + fundamental doublet Higgs theory 

L = \f;^F;^ + {D,<t>y{D,<t>) - ^m'<t>^<t> + A(</>t</,)^ (14) 

We have thus omitted from the full standard model the U(l) subgroup as well as 
fermions, which - with some reservation for the top quark - are inessential for the 
thermal phenomena studied. In fact, fermions can be trivially included in the scheme 
of dimensional reduction. Their effect is to modify expressions relating 3d parameters 
with 4d ones, but they do not in any way change the qualitative nature of the effective 
3d theory. 

The leading principle of our analysis is the replacement of the full 4d theory (14) by an 
effective 3d theory, obtained by integrating over the non-static modes perturbatively. 
As we discussed in the previous section, this is a very good approximation to the 
problem under consideration. 

Integrating over the non-static modes to 1-loop accuracy in the MS scheme the 
following effective action is obtained: 

5eff[4"(x), AS(x),</>i(x)] = I d'x{^FtjFtj + ^{DiAoTiDiAoT + {D4)\D4) + 

+ \mlAlAl + \\A{AlAlf + mlct>^ + Uct>^f + h^AlAlct^^^ (15) 

The key property of the theory (15) is that all the couplings Xa, A3, and gl are 
3d renormalization group invariants (there is no ultraviolet renormalization of them). 
These couplings are given in terms of the 4d couplings by 



As = T 

1 



gl = g\fiT)T, (16) 
K^^T) + rr^.lg'M], (17) 



4^3 



167r2 8 

1 / 49 1 M 

1 + (12A(/xt) + y/(/^T) - :^9\fiT))\ , (18) 



A. = ^^T^, (19) 



17g\fiT)T 
487r2 

where the 4d couplings are run to the scale fiT by the standard f3 functions. 

On the contrary, the 3d mass operators for the Higgs and Aq fields contain linear and 
logarithmic divergences. To renormalize the theory, one adds to eq.(15) the counterterm 
action 

SS=- j + ^^ml^A^Asj, (20) 

where 

Srm? = /imSlin + /2mSlog, Sm'^jj = /iLlSlin + /2LlSlog. (21) 
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Here 



and 



/im = ^^3 + 6A3, /lL> = 5(5-3 + A^), 

81 

f2m = -gt + 9Ws - 12A^, f2D = 5X1 - 20glXA 
Id 



(22) 



(23) 



are exact expressions with no higher-order corrections and and Siog are counter- 
terms depending on the regularization scheme used. The linearly divergent counterterm 
Siin vanishes in dimensional regularization; on the lattice it behaves like = S/47ra 
(S = 3.1759114, a is the lattice spacing). The quantity Siog is related to the log- 
arithmically divergent 2-loop sunset diagram (eq.(22) of [2]). In the MS scheme it 
is 



and on the lattice 



^log 



OC 



IGtt' 



647r2e 



log(ar). 



In the MS scheme the masses are scale-dependent 



which can be integrated to give 



1 dmKus) 

J 2m 5 /^S " 



IGtt' 



2D 5 



"^sC/^s) = T^/2mlog— , m^/is) = 2 
lOTT'' fls lOTT'' 



1 . ^ Ad 

f2D log , 



(24) 
(25) 

(26) 
(27) 



where and Ad are integration constants. The constant can be obtained by 
comparing the 4d and 3d calculations: 



^^^^ + 2^^^^^ + W.^ V 96 



167 , 3, 



-m 



2 ' ' 167r2 

gl A49 2 , 3 



3r 

/2m (log he 



(28) 



1 2 1 
—m„ H - 

2 ^ 167r2 



/2m ( log — + c 



where m|j = m'^[fiT) and [26] 

c = [log(87r/9) - 27£; + C'(2)/C(2)]/2 = -0.348725. 

For rriD one can similarly write 

5 1 1 

"^i)(/^3) = -g^{nT)T^ + 7^/21? log( — ) + -r^icigt + c2glX3), 



(29) 



(30) 



but the constants Ci,C2 and, hence. Ad have so far not been computed. 
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In practice, one can write 
2 



As 



8m^ ' 



mw^ = 80.6GeV ha = -gl. (31) 



The explicit form of the total 2-loop effective potential is needed frequently in what 
follows. The tree and 1-loop parts are 



Vi{<t>) 



(32) 
(33) 



and the 2-loop part is (eq.(33) of [2]) 



35-1 ,2 



2iZ"(mi,mr,mr) - -iZ"(mi, m^, 0) + iZ"(mi, m^, m^) 



m 



I 5-[-H"(mi,mr,0) - iZ"(mi, m^, m^)] 



m 



4m^ 



[if(mi, 0, 0) + H{rni,mTi rriT) — 2H{mi,mTi 0)] 



m-. 



mi 



2mT A.m\ 



"iX^cj) [iZ"(mi,mi,mi) + iZ"(mi,m2,m2)] 



63 -, 3 - , ,41 

—H{rnT,mT,mT) + — iZ"(mr,0,0) - — 



zgllif^T ~ 4:1711) H{'mL,mL,mT) ~ ^mTrriL - m\] 



+4:glm^ + -gl{2mT + mL){mi + 3m2) 
o 

15 3 
+— A^m| + -Xaiml + 2mim2 + Smj) 



8 



fi'siC^^r ~ — 2ml)H[mi,m2,mT) + (wiy — Ami) H {^1712, m2,mT) 



{ml - mlf 



mn 



■[H{mi,m2,mT) — H{mi,m2,0)] 



+(m^ — ml)[mi — m2)/mr + mr(mi + 3m2) — mim2 — "^2]|- (34) 



Here the masses are 



■mr = -g3<t>, 



2 2 I -'- 2 /2 



"^1 = "^sC/^s) + 3A3^^, ml = ml{ii2,) + Xzcj)^ 



(35) 
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and the finite part of the sunset function is 



H{mi,m2,m3) = log , , h ^. (36) 

The coefficient of |log(/i3)^^, coming entirely from the if-terms, gives the value of 
/2™/167r2. 

This potential has been calculated in the Landau gauge. In a general covariant 
gauge the potential has been computed in [27, 28]. The 4d computation of the high 
temperature asymptotics of the potential [20]-[22] leads to the same result as the 3d 
one [2]. 

The 3d scale fis appearing in (28) and in (34) is arbitrary. We refer here to [2] for a 
discussion of the renormalization group improved effective potential, which allows one 
to fix this scale in a natural way. 

4 Gauge-invariant order parameters — condensates 

As was argued in ref. [2] a reliable perturbative computation of any physical observable 
in 3d theory should be done at least on the 2-loop level. The reason for this is the fact 
that the logarithmic renormalization of masses in 3d starts at the 2-loop level. The 
only quantity known at present for a 3d gauge- Higgs system on 2-loop approximation is 
the effective potential in Landau gauge [2] or in arbitrary gauge [27, 28]. The effective 
potential, is, however, a gauge-dependent quantity, and cannot be immediately used 
for the extraction of gauge-invariant information. In this section, we show how simple 
manipulations with the effective potential can generate a set of gauge-invariant order 
parameters - condensates. Those condensates can be computed with 2-loop accuracy 
with little extra work and can be easily measured on the lattice. 

By condensates, in general, we mean vacuum averages of any local composite gauge 
invariant operators, 

_ /^^e-^WQ(^) 

Examples are provided by the scalar condensate {<f)^<f)), the "gluon" condensate [F^-F^-), 
etc. These composite operators contain the product of fields at the same point. There- 
fore, their expectation values are divergent (infinities cannot be removed by countert- 
erms renormalizing the masses, couplings and fields of the underlying theory). So, their 
value is prescription dependent and cannot be fixed unambiguously. 

We would like to provide a natural prescription for a number of condensates related 
to the ground state energy. We define its renormalization prescription in the next 
subsection. 

4.1 Renormalization of vacuum energy 

The renormalization of our 3d theory requires the introduction of 2-loop mass countert- 
erms for and m|,. After that, the computation of any physical amplitude does not 
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contain ultraviolet divergences. However, the computation of the value of the ground 
state energy defined by the functional integral 



y"x»^e-^('^), (38) 



with V3 being a spatial volume, does contain ultraviolet divergences. We extend the 
prescription of the MS scheme (remove poles in e) to the vacuum graphs. Thus, vacuum 
counter-terms are 

lD7r"'4e 

^4'* = -Jl^^^^st + 6253^3 + 63^3^1 + e4A^], (39) 

(47rj*4e 

(Aj4-terms are, for brevity, not included). Since in the MS scheme linearly divergent 
integrals are equal to zero, there are no 1- or 3-loop counterterms. After this renor- 
malization the vacuum energy density e^ac is finite, but fis dependent. For example, 
the 2-loop expression for the ground state energy of the unbroken phase found from 
Vms{(I> = 0) in (34) is: 

V{0) = hVr{0) + hW2{0) 

^ -[3ml + Amlifis)] 



l"/r» 2 1 ^2 ^ \ \ 2 2 



167r 



2 



(40) 

Clearly, the fis dependence here is a reflexion of the arbitrariness in the definition of 

The relation of the ground state energy to the effective potential is obvious. The 
latter is defined through 



^-V3V{4>)-J4> 



|x)^e-^W-^'> = e-^(-^), (41) 



with (j) = W'{J). The external current J is needed to give the scalar field the required 
value (j). If the external current vanishes, V'{(j)) = —J = 0, the system settles in its 
physical ground state and by evaluating eq.(41) we obtain the value of the effective 
potential at its minimum, V{v[T)) , V'{v[T)) = 0. This is precisely the ground state 
energy Evac of the 3d theory defined above by (38). 



4.2 Renormalized condensates in MS scheme 

Now we are ready to define gauge-invariant condensates. Our 3d theory contains five 
different parameters, , mg , m|, , A3 and Xa- So, there are five special condensates 
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that can be computed by a simple differentiation of the ground state energy density. 
Since Evac is finite, all its derivatives are finite (but /is-dependent). So, we define the 
renormalized condensates as corresponding derivatives of e^ac- They are related to the 
unrenormalized condensates as follows: 
1. Quadratic scalar condensates: 



2 omjj 2 D47r 



e 



Quartic scalar condensates: 



3. Gauge condensate: 



4l^i.-^i.;ii-53 ^^2 - 4\^i,-fi,;+53 ^^2 64^2, 



All renormalized condensates defined above are finite but /is-dependent. 

For completeness, we also write here the relations between the condensate of the 
kinetic part of the action with other unrenormalized condensates. These follow from 
the independence of the ground state energy on the normalization of the (j) and Aq 
fields in the functional integral (38): 

{{Di(t>)\Di(t>) + {ml + 8m'')(t>H + 2\^{(t>Hf + h^AlAl(t>H) = const, 
{{DiAl)\DiAl) + {ml + 6m')AlAl + MA^A^f + 2/i3AS^S</>V) = const,(45) 

which are nothing but Schwinger-Dyson equations. The constants here do not depend 
on the parameters of the theory. 

The condensates themselves do not have much physical meaning, just because they 
(and the ground state energy) are dependent on the normalization point. However, 
if the system exhibits a first-order phase transition, then in the differences between 
condensates in different phases the fis dependence cancels out. An important physi- 
cal characteristic of the system, namely the latent heat, is related to the jump of the 
renormalized (or, what is the same, unrenormalized) Higgs scalar condensate. Another 
application of condensates is that they can be measured on the lattice with high accu- 
racy. Below we will provide a relation between condensates on the lattice and in the 
continuum. 
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4.3 The latent heat of the phase transition 

The computation of the latent heat of the first-order phase transition is based on three 
equations. 

1. The ground state energy density depends on six dimensionful variables and has the 
dimensionality GeV^. From here we get 



+ + >^A^7- + 2ml-— ^ + 2"^I)^-^ + /^s^ — ]£vac = 3e„ac- (46) 

Ogi 0X3 OXa Om^ Omjj O/ia 

2. The difference between the ground state energies Ae of the broken and unbroken 
phases is renormalization group invariant (/is-independent). Therefore 

dAe^ac hrn f^D 

dfia Ofia IGtt^ Om^ IGtt^ Omjj 

3. The derivative of Ae with respect to the temperature for g,X and tuh fixed at 
the transition point (where Ae = 0) is just 

r^^ = <rA(<^V). (48) 

Remembering that, in 4d notation, Aevac = Ai^/T = —Ap[T)/T, this is nothing but 
the latent heat L = A[Tp'[T)] of the transition. Here F is the free energy and p is 
the pressure of the system. According to (42), the jump of this order parameter does 
not contain any divergences, and, therefore, does not depend on the regularization 
scheme. The simplicity of expression (48) could be contrasted with the form of the 
corresponding expression in 4d theory [6]. 

Another quantity which may be of interest is the difference between the ground state 
energies in the close vicinity of the phase transition. It is given by 

Ae_ = ^ [\Fr,Fr, - 2{ml + Sm^ + ^)4>^ - X,{4>'4>f 



-{ml + 8ml + - -M^l^lf 



). (49) 



We stress that all condensates appearing in this expression are unrenormalized, so that 
there is no difficulty in putting it on the lattice. 



5 Perturbative results for the ground state energy 
and 

In this section, we shall compute the ground state energy e^ac = ^{"^iT)), V'{y{T)) = 
and the scalar condensate (^V) in 2-loop perturbation theory. The second can be 

^^We neglect here the slow logarithmic variation of all coupling constants with the temperature. 
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obtained from the first by the use of eq.(42) and the first can be obtained from the 
2-loop effective potential (34). 

In perturbation theory the effective potential is derived as an expansion with respect 
to the Planck constant which serves as a loop counting parameter: 



N 



y(</>) = 5:rK(</>). (50) 



n=0 



Simultaneously, the equation for the determination of the vacuum expectation value 
v[T) is 



= 0. (51) 

viT) 



n=0 d(f) 

We should now solve v[T) from (51) and insert in (50). Two different regimes should 
be distinguished. 

1. The "classical" regime. In this case, spontaneous symmetry breaking occurs on 
the tree level, and the equation dVo{(j)) / d(j) = has a non-trivial solution, v'^[T) = — 



If the expansion parameter p = 5'|/(7rmr) = gl^3 / {—fnl) is small enough, eq. (51) 
can be solved perturbatively: 



v{T) = j:h-v^{T). (52) 

n=0 

Inserting this in (50), the ground state energy becomes a power series in h up to the 
order N: 

N 

£vac = 'Y K^Vn{<f))\4>=v{T), (53) 
n=0 

where all terms containing powers of h higher than N should be dropped. It is assumed 
that Vn{(j)) are re-expanded in powers of h. This approximation works provided that 

ml < 0, \ml\ > glXs. (54) 

2. The Coleman- Weinberg regime. Here spontaneous symmetry breaking occurs 
due to radiative corrections and the equation dVo{(j)) / d(j) = need not have non-trivial 
solutions. The conditions (54) are not satisfied. In this case a perturbative solution of 
eq. (51) makes no sense. So, eq. (51) should be solved exactly. We denote the solution 
by vn[T). Again, if the parameter p = gl/[TrmT) = 2g3/[TrvN[T)) is small enough, the 
quantity 

N 

e.ac = E^"^4^iv(r)) (55) 

n=0 

is a good approximation to the ground state energy. One usually enters the Coleman- 
Weinberg regime when the tree scalar mass is small enough and when A3 <C gl- For 
ml = 0, expression (55) provides an expansion of the vacuum energy with respect 
to Xa/gl; in each order of Xa/gl the summation of all powers of h is automatically 



16 



performed. Very close to the phase transition the first method certainly fails and the 
second method should be used. On the contrary, if condition (54) is satisfied, then a 
Coleman- Weinberg type of computation of the ground state energy can be performed 
as well. In 3d, if the effective potential is computed up to the terms of order h , then 
the difference between the two computations is of the order of h^^~^^\ The fractional 
power of h is due to the contribution of Goldstone bosons to the effective potential, 
which produces infrared- dangerous terms in the second type of computation. Since 
fractional powers of h must be absent, the first method of computation should be used 
in this regime. 

Below we shall derive explicit formulas for the ground state energy in the broken 
phase in the 2-loop approximation in the classical regime. The Coleman- Weinberg type 
of computation does not require any additional analytic work: simply, one numerically 
minimizes the 2-loop potential defined in eq.(34) and then determines the ground state 
energy at the minimum of the potential. Finally, one can use (42) for the determination 
of condensates. 



5.1 The ground state energy to order h 

Let us compute Evac = V{v[T)) in the classical regime, when spontaneous symmetry 
breaking occurs already at tree level, 

= (56) 
It is assumed that —ml > 0, and sufficiently large. At this point the masses are 

- 2 ~"^35'3 -2 2,-2 -2 o2 - 2 n 

my = , m^ = mjj + m^, m-^ = —Zm^, mj = U (57j 

4A3 

and the leading approximation to V{v) is: 



V.(v) = -g. (58) 

The corrections to eqs.(56) and (58) can be obtained from the 2-loop result (1.33) 
for the effective potential as follows. The definition of v is 

V'{v) = 0, V{v) = Vo{v) + hVr{v) + h''V2{v) + 0{h^), (59) 

where v"^ coincides with only to leading order. To higher orders (note that v^^^ is 
the order contribution to v^, not (vnY) 

= i;fo) + + ^'4) + 0{h^). (60) 

Inserting this in V^v"^) = (prime means derivative with respect to ^^) and expanding 
gives the equations 

"^3 + M'o) = 0' (61) 
M'i) + 2K(4)) = 0, (62) 

^3^) + 2KV(o)H) + 2^2H)) = 0- (63) 
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For later use, we tabulate the partial derivatives of Vi'. 



dVi 
dml 
dVi 

d(j)^dml 
From these one finds 



h , 

— (mi + 31712), 

OTT 

3^ \ r 33 / c\ \i 

A3["^l + "^2 + 7T-(2mr + mL)\, 



Stt 
16^' 



■A, 



1 1 A 

— + — , 
mi 1712 J 

3 1 

— + — + 
mi m2 



9l 



4A, 



- + -) 

mr mi/ 



^(0) 



'(2) 



^^0 



-m 



A., 



A 



3 

47r 



53^ 



mi + — — (2mr + m^) 
4A3 



(64) 
(65) 
(66) 
(67) 

(68) 
(69) 



327r2 



3 1 

— + — + 
Lmi m2 



9l 



4Af 



2 1 

— + — 
mr mi 



["^1 + 'T^C^'^T + mi)] 
4A3 



|-^2Wo))- 
^3 



(70) 



We have here explicitly written out the "disconnected" parts that follow from the 
derivatives of the 1-loop potential. The 2-loop expressions follow automatically from 
eq.(34), but are much lengthier. Note that the second derivatives of Vi in (66,67) 
contain terms ~ l/m2, which diverge at the saddle point. Similar terms appear in the 
first derivatives of V2. Evaluating the first derivative dV2/d(j)'^ one can show that the 
l/m2-terms cancel between the 1- and 2-loop parts in eq.(70) before taking the limit 
m2 — > 0. This cancellation actually takes place only in the Landau gauge [27], but this 
is not surprising since v[T) is not a physical gauge invariant quantity. 
With this expansion of v'^[T) we can write 



V{v{T)) 



yo(4) + hvl^ + E'^)) + hv^ivl^^ + hvl^) + ^2^2(4)) 



4A3 



127r 



(6my + 3m|^ + mj) 



(71) 



V2("^^,"^L,"^l,"^2 = 0) 



9A3 
647r2 



mi + -^(2mr + m^) 
4A3 



In other words, the value of V{v) in the loop expansion is given by the value of the 
2-loop potential for the saddle point masses in eq.(57) corrected to order h"^ by the last 
term in eq.(71). This corresponds to a set of 1-particle reducible diagrams, which are 
neglected in the computation of the effective potential. 
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We shall later also need the 3- and 4-loop contributions to V{v[T)). These are 



(72) 



where all the potentials should be evaluated at the saddle point as in (71). These 
equations can be used for an estimate of the higher-order corrections to the effective 
potential through the lattice measurements of the scalar condensate. 

For completeness we also present 2-loop expressions for the unbroken phase. If 
ml > 0, there is another symmetric phase saddle point at = with the associated 
mass values 

niT = 0, niL = m£), m\ = ml = ml (73) 
and a vanishing leading minimum value 



Vo{0) = 0. 



(74) 



At the point (73) higher-order loop expansion will lead to infrared divergences, but 
to this order one similarly obtains v = 0[h^) = ((^V) is non- vanishing, see below) 
and V{0) = hVi{0) + E^V2(0) as given by eq.(40). Note the explicit /is-dependence 
here. As it should, the difference V{v[T)) — V{0) computed from eqs.(71) and (40) is 
independent of fis to this order in h. 

5.2 Computation of (^V) 2 loops 

The result for {(f)^(f)) to 2 loops can now be directly obtained from the expansion (68-70) 
and the general formula (42). According to (42) we have, in the broken phase (/ (Px/Va 
is implied in the averages). 



(</>V) 



5(^0 + ^1 + ^2) 



dml 



d 



+ 



(0)^ 

dmld(j)^ 



5^2(4)) 



m% 



(75) 



where 



(</>V)(i) 



dml 



(</>V)(o) 



-m% 



2A3 

1 



"^1 + + "^i) 

0A3 



(76) 
(77) 



(</>V)(2) 



2^2) + 



dmldcf)^ 



4) + 



dml 
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1 d'^Vr ^ g'Fi \ / ]^dVi\ ^ ( dV2 dV2 



2 



A3 5</>25</>2 dmld(l)y\ Xsd(l)y \dml X^dcf) 



["^1 + 'T^i'^'^T + mi)] 



All derivatives here should be evaluated at the saddle point masses in eq.(57); in 
practice those of V2 are most conveniently evaluated numerically. The 1-loop expression 
(77) will be written later in a more general form, appropriate for lattice regularization, 
by using the tadpole function I{m) introduced in Appendix A. Note that here also poles 
in m2 are generated by the derivatives dm2/dml = l/2m2 and dm2/d(j)'^ = A3/2m2. 
Because ml = ml + As^^, these precisely cancel in eq.(78) and leave only the terms in 
the last of eqs.(78). As shown in [27] this cancellation takes place for any gauge fixing 
parameter ^, as should for a physical quantity. 

In the symmetric phase one similarly obtains from eq.(40) that^^ 

(<^V)s._ = + 1^3 log ^ - + 2A3 + . (79) 

In a 3d theory the perturbative computation of {(f)^(f)) contains a linear 1-loop and 
logarithmic 2-loop divergence and these have been treated above in the continuum MS 
scheme. Due to their importance for the lattice computation it may be useful to specify 
their origin once more. Their coefficients can be obtained as follows: 

(2<^<^) = 2 ^-^ + (E'^.)div 

dV{v{T),ml) , r d'p 1 



+ 4 



/ 



dml J (27r)3 + Ap 

JV{v{nml) f d'p 1 3gl f d'p 1 

where A = —3gl/W is the coefficient of the linear term in the 1-loop Higgs self-energy 
(n^(A;) = Ak, [2], eq.(83)). In other words, the linear divergence comes from the simple 
scalar loop and the logarithmic divergence from the linear term in the scalar self energy 
caused by the emission and absorption of a gauge particle. 

The computation of the vacuum expectation value of the composite gauge invariant 
operator (f)^(f) is related to the effective potential for this quantity introduced in [13]. 
We clarify this relation and compute the corresponding effective potential to 2 loops 
in Appendix B. 



^^We would like to note here that since perturbation theory does not work in the unbroken phase 
this expression cannot be used for any comparison with the lattice results. 
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5.3 Computation of {AqAq) to 2 loops 

In contrast to the doublet Higgs field ^, the triplet scalar field Aq is dynamically 
inessential and can be integrated over [2] . Due to the relatively large value of the mass 
m£), the vev of the triplet scalar field is zero. For pure SU(2) gauge theory this has 
been numerically studied in [29]. Thus the field feels the phase transition only in 
that the field (/) in m| = m|, + g'^(/)'^ /A develops an expectation value v[T). In analogy 
with (j) in the symmetric phase (eq.(79)) in 2-loop perturbation theory one can write 
for the finite part 



dV{v{T)^l) _ 3 



3 f 2 



16 



/is 2mij — 2mi — mTrriL 
4 log h 



SrriT + mi + 3m2] 5 



rriL 

The linear and logarithmic divergences are 



}• (81) 



KA;U,=3/p^- + f /^-. (82) 

6 The lattice action and the curves of constant 
physics 

6.1 The lattice action 

The lattice action corresponding to the continuum Lagrangian in eq.(15) has the fol- 
lowing form: 

^ = /5GEE(l-^TrP,,) + 

+ ^/5GEE[Tr^o(x)C^-^(x)Ao(x + i)C/,(x) - Tr A^(x)] + 

+ E/52^^Tr A^(x) + E/5f (^Tr A^(x))^ + (83) 

E Et^Tr $t(x)$(x) - ^Tr $t(x)C7,(x)$(x + i)] + 
+ YXi^ - 2/5« - 3/3i^)^Tr $t(x)$(x) + /^^(^Tr $t(x)$(x))^] 

X 

-^/5^E[^TrA^(x)lTr$t(x)$(x)]. 
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Here Ui{'x.) and Pij are the standard link and plaquette variables, the scalar field $ on 
the lattice is related to the continuum scalar field (j) through 

^ = VRl, i2i = |^</>V = ^Tr$t$, (84) 

where a is the lattice spacing, ^ is a unitary SU(2) matrix, Rl is the radial mode of 
the Higgs field, and the lattice matrix field Aq is given in terms of continuum field Aq 
as ^ 

^0 = -igsaTaA^. (85) 

The main aim of this section is to find the relation between the parameters fSn, 
(3r] (3^ and (3^ of the lattice action and the parameters A3, m^] Xa-, m\) of the 
3d continuum theory. The super-renormalizable character of the 3d theory allows one 
to write the relation in a form that is exact in the continuum limit iV — )• 00, a — > 

0. N I fia 00. The 3d continuum parameters were related to the physical 4d ones in 
the previous section. 

The lattice-continuum relation is established in several steps: 

1. In subsection 6.2 we formulate the general problem of defining the constant physics 
curves in the 3d gauge- Higgs system. Then (subsection 6.3) we consider simplifications 
arising in electroweak theory and show that the problem is reduced to the computa- 
tion of three pure numbers [rj^fj^rj) - originating from three different classes of 2-loop 
diagrams contributing to the Higgs self-energy. In subsections 6.4-6.6 we compute two 
of those numbers [fj and rj) analytically and formulate a way of MC computation of 
the third one (?/). 

2. In subsection 6.7 we work out an exact relationship between the condensates in the 
lattice regularization scheme and MS scheme. 

4. Finally, in Section 7 we use the 2-loop continuum computation in Section 5 and the 
relationship between different renormalization schemes in subsection 6.7 to determine 
the constant rj by comparison of the lattice data with perturbation theory. This allows 
one also to estimate the magnitude of 3-loop corrections to the effective potential and 
the magnitude of 2-loop finite a effects. 



6.2 Constant physics curves: general formulation of the prob- 
lem 

In the tree approximation all the five lattice coupling constants are given in terms of 
continuum ones by the following equations, which directly follow from the discretization 
procedure and the form of the continuum Lagrangian: 

4 1 

f^G = --, 86 

f3n = \\saf3'H = (87) 

f^t = H^o, (88) 

93 
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_ 2(1 - 2Pn - .89^ 

m3|tree- , (89) 



2 I r, f^2 



(90) 



The first three relations (86,87,88) are not spoiled by the renormalization procedure 
and are exact in the continuum limit, while (89) and (90) are modified when radiative 
corrections are taken into account. As we have discussed in [2], 1-loop and 2-loop 
counterterms must be added, so that in the continuum limit a — )• we have the exact 
relations 

— ^ = m^i^is) - hm- — ^ log( ) - -— , (91) 

Pho, 47ra Idtt"^ a/is Idtt"^ 

-2-5-^ = mJJ{^l3) - fiD- TT-^ log( ) - -— , (92) 

PgO' 47ra Idtt"' a/is Idtt"' 

defining the constant physics curves. Here S = 3.17159 (Appendix A) and the coef- 
ficients fi of the linear 1-loop and logarithmic 2-loop terms were given in the section 
3. However, one also needs the 2-loop constant terms (when a — )• 0) and P2D 
and their determination is our main task. Note that the above relations are, in fact, 
/ts-independent, due to the /ts-dependence of the mass parameters in the minimal sub- 
traction scheme. 

Exactly as /2m and /2D, the constants P2m and P2D are, for dimensional reasons, 
quadratic polynomials constructed from the coupling constants 5'3,As and Xa- The 
coefficients in the polynomial Pi can be determined by explicit computation of the 2- 
loop mass renormalization in lattice regularization. There are at most six dimensionless 
numbers (rji) determining P2m'- 

P2m = Vgdt + ^aAs + TJaXa + VgXal^?. + VgAgl^A + VAX^s^a- (93) 

From the structure of the 2-loop scalar mass renormalization graphs, one can see that 
f]A\ = f]A = ^1 so that four numbers should be determined. In the Aq case we have, in 
complete analogy: 

P2D = Pggt + P\^l + PaX\ + PgXgl^S + PgAgl^A + PAxXsXa- (94) 

For the same reasons as in the scalar doublet case we have pA\ = p\ = 0. So, the 
general analysis of this 3d theory requires a knowledge of eight pure numbers which 
fix the relation between the lattice renormalization and MS schemes. To find them, 
one should compute the renormalization of the Higgs and Aq masses on the 2-loop 
level in MS and the lattice regularization schemes, and compare the results. Once they 
are known, there is a one-to-one correspondence between the lattice couplings and 
continuum parameters which allows one to relate the results of lattice simulations to 
the 3d effective theory. Then, to relate 3d and 4d high- temperature physics one should 
use (28) and (30). The two so far unknown constants in (30) are to be computed by 
comparison of continuum 3d 2-loop renormalization of the Debye screening mass with 
corresponding 4d computations. To summarize, we have 10 numbers to be determined 
to relate 3d lattice simulations to 4d high temperature theory. 
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6.3 Simplifications for EW theory 



In the context of the EW transition the problem can be simplified, with a good accu- 
racy, by making use of the fact that the coupling constant Xa = ^^dl is numerically 
very small and that there is no ^f^-term in f2D- This allows one to use a formal 
expansion with respect to the 4d coupling constant assuming for power counting 
that A ~ and keeping at most the terms ~ g'^ in the previous expressions. In 
this approximation^^/2D = 0, while the polynomials and P2D are defined by five 
constants ?7g,?7gA5?7A and Pg,Pg\- In more detail, we shall parametrize 

81 

P2m = —gtv + 9A353^ - 12A^i7 + f2mC, (95) 

where c is the constant defined in eq.(29). This will cancel when mKfis) is related to 
4d physics. 

Another simplification takes place when we notice that the mass of the Aq field is 
parametrically and numerically larger than that of the vector boson and scalar field 
near the phase transition in the broken phase, at least for a sufficiently heavy Higgs 
boson [2]. In addition, the Aq mass in the unbroken phase is parametrically larger than 
the typical 3d mass scale gl- So, in principle the Aq field can be integrated out and 
a theory without the Aq field with modified coupling constants can be formulated [2]. 
From eqs. (1. 51,1. 52,1. 53) one can see that the uncertainty of the order of g'^ in the Debye 
mass gives an order g^ uncertainty in the gauge and scalar self-coupling constants and 
an 0{g^) uncertainty in the Higgs mass. These are higher-order terms according to our 
convention. Therefore, we take P2D = and omit the 0{g'^) terms from eq.(30). 

Note that the formal suppression of the higher-order terms does not necessarily imply 
numerical suppression. We consider the systematic uncertainties associated with the 
ignorance of those terms below. 

To summarize, the approximate constant physics curves for the electroweak theory 
are given by eqs. (91, 92) with 

f2D = P2D =Va= VgA = VAX = (96) 

and P2m given by eq.(95). 

In spite of all these simplifications we still have to compute three numbers, t/,?/ and 
rj. In the next subsection and Appendix A we carry out an analytic computation of the 
parameters rj and fj. The one remaining parameter will be determined by a combination 
of analytical and Monte Carlo methods. 

The previous discussion referred entirely to the 3d effective theory. To have an 
explicit relation between the temperature and lattice parameters, we introduce into 
eq.(91) the explicit relations between the parameters of the 4d and 3d theories in 

^■^Formally, f2D is multiplied by log ^ = log ^ , which is singular in the continuum limit. How- 
ever, in real lattice simulations with Pa < 40, the value of this log never exceeds 2, and this term may 
be neglected. We will estimate the possible influence of this term on the results below. 
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eqs.(28). Using the parametrization (95) the constant physics curve becomes 
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It is illuminating to define from here fSn for large f3Q. One finds that 

2mH ^2 



MT) = Ph{T = 00) + {j^) + 0(/35'), 
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7--(A-p/12)(A-p/6)-5 +0(/35'), (101) 



where the abbreviations are: p = m'^/m^, A = 9S(1 + p/3)/8ir and B = second 
line of eq.(99) without sign and the factor l/jS^. ne sees that the T-dependence of 
fSff is very simple and that the interval between fSn = 1/3 and fSn = f^niT = 00) is 
unphysical. 

The constants ?/, fj and rj can be found by computing the 2-loop Higgs mass renor- 
malization in the lattice scheme. Instead of considering only this, we prefer to work in 
an equivalent language, namely with the effective potential. This will allow us to get 
simultaneously a number of results concerning the values of different condensates in 
the ground state. To make the method clear we shall first consider the simplest case, 
the computation of rj. 

For those who are not interested in details of computations we give the result below: 



77 = 2.18(6), 77 = 1.01,77 = 0.44, 
where only MC statistical errors in determination of 77 are shown. 



(102) 
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In [4] some combinations of these numbers, relevant for the study of niH = 80 
GeV Higgs boson were presented. At that time the quantity f] was computed by MC 
methods, rather than analytically. Due to a mistake in the computation of parameter 
7/, made in [4], the numbers in [4] must be corrected as follows: in eq.(ll) of [4] rj = 2.0 
instead of 1.35, and in eq.(13) fj = 0.66 instead of —4.70. This change slightly affects 
the extraction of the critical temperature from the lattice data, but does not affect 
qualitative results and conclusions of [4], see [3]. 

6.4 The 2-loop lattice effective potential in scalar theory: 
computation of 

The constant fj multiplies the scalar coupling constant. Therefore, it can be computed 
in pure scalar theory with SU(2) symmetry. To be more general we take the SO(N) 
scalar theory with the scalar field in a fundamental representation, considered in Ap- 
pendix B.l of ref. [2] (the case we are interested in corresponds to iV = 4). Since the 
perturbative computations in scalar field theory on the lattice are not so complicated 
as in gauge theories, the 2-loop lattice effective potential accounting for finite size and 
finite a effects can be quite easily computed. 

The continuum Lagrangian is (in this subsection, A3 — )• A,m3 — )• m) 

jc = l{d.<t>ay + lm'<t>i + \\{<t>iy, (103) 

and the 2-loop effective potential for the scalar field ^0 in the MS scheme is: 

Vms(^^o) = ^m'^{n)(f)l + ^\(f)l - 

-7^[ml + {N-l)ml]+ (104) 
+ + 2(iV - l)mim2 + {N^ - l)m^] - 

_3 J!- flog ^ + log " 

167r^ V 3mi 3 mi + 2m2 6 / 

where ml = m'^[fi) + 3A^0 5 "^2 — "^^(/^) + '^^^o 

/^^^ = -t^, hm = -2{N + 2)X\ (105) 
Ofi 167r"^ 

In order to find the relation between the MS and the lattice regularization schemes, 
we have to perform the corresponding computation on the lattice. We take the lattice 
Lagrangian in the form 

= ^{A,<t>.y + ^myi + \\{<t>iy, (io6) 
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where is a lattice difference in the z-th direction and ms is the bare mass. Then 
the effective potential is given by [h is a loop counting parameter) 



wi 



th 



i^l-loop = J{miL) + {N- l)J(m2L), 



(107) 
(108) 



i^2-ioop = ^mrmLY + 2{N-l)I{rmL)I{m2L) + {N'-l)I{m2LY] 

N — 1 

-3A^[ifL(miL,miL,miL) H — ^fL("^lL, "^2L, "^2L)]^^05 (109) 

where m^^ = + 3A^q, = + A^q, and the functions 7, J and Hl are the 
lattice analogues of the continuum tadpole, 1-loop energy and sunset diagram functions, 
defined in Appendix A. For small am these functions have the following expansion in 
powers of the lattice spacing a: 



I{m) = [S - (am) - ({amf + 0{{amf)\, 



(110) 



= 7^[const + ^S(am)' - l{amf - \^{amY + 0{{am)% (111) 
47ra Z o 4 
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where the numbers ^ and ^ were determined numerically, ^ = 0.15281, C, = 0.09. 
Expanding for small a one obtains 
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With the choice of the following relation between the lattice and continuum masses 



ml = m'U) - hX(N + 2)^ + h'^'^^^ 
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we get 



Note, in particular, that 



NT, 

Vl = Vms + h^m\ 



(115) 
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+ 0{h^) = ml + 0{h^), (116) 



that in 2-loop terms rriiL and mi are equivalent and also that the term ~ am\j^ — > 
constant when a — )• 0. Of course, is //-independent. In terms of the parametrization 
in eq.(95) we have 

Plm = f2m{V + C), fj = C-C = 0.44. (117) 



This completes an estimate of fj. 



6.5 The effective potential in MS and lattice regularization 
schemes 

Since the lattice regularization provides nothing but another subtraction scheme in 
perturbation theory, and since the only divergences of the 3d theory are those related 
to mass and vacuum energy renormalization, the exact lattice effective potential can 
differ from the exact MS continuum one by only two types of terms. The first one, 
multiplying ^^/2, is associated with the Higgs mass renormalization, and knowing 
it would give us the unknown constants in P2m- The second is a ^-independent piece 
connected with the renormalization of the ground state energy. Just by power counting, 
it comes from diagrams up to 4-loop order. It is irrelevant for the study of the phase 
transitions but will be related to the different condensates introduced in Section 4. 
Generalizing the result for the scalar theory derived above, the a — )• limit of the 
lattice effective potential is thus 

VL{<j>) = Fms(</>) + /i™^^</>' 

+ T^[/2™(log — + c) + ^gtv + 9\sgifj - 12\lij]U' 
loTT'' afi Id 2 

+Vl^\ (118) 

In section 4 we have already specified the renormalization of the vacuum energy or 
the convention for the value of Vusifl^ = 0)- This is a nonphysical divergent quantity 
and its value could thus be fixed at will. Now when it is fixed, any additional lattice 
effect is contained in V^^*^. 

The general form of V£^'^ to different orders in the loop expansion can be fixed as 
follows. The 1-loop contribution to it, V"^™*^, can be found from the 1-loop lattice 
effective potential, which as a direct generalization of eq.(108) or of eq.(33) is 

Vilop = 6J(mr) + 3J(mi) + J(mi) + 3J(m2), (119) 
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where the particle masses are given by eqs.(35). To this order the lattice bare mass 
and ml are equivalent. Including the m^/a-term in the expansion of J{m) one has 

V^r = ^ml + f^ml. (120) 

The 1-loop term comes just from a bare loop and no coupling constants enter. To 2 
loops, first powers of g\ and A3 enter and the structure of the 2-loop constant term is 

= T^,[rnl{a^gl + a^As) + ml{a^gl + a4A3)]. (121) 

From the consideration of possible 2-loop diagrams, 02 = 0. Similarly, from the compu- 
tation of the previous section we find directly that 04 = (in scalar theory the 2-loop 
lattice potential for a — )• vanishes when ^ = 0, see (115)). The dependence of 
and as on the renormalization scale can be established from the obvious fact that the 
lattice effective potential or VM's{(t> = 0) + V^^"^ must be /is-independent. From the 
explicit expression for VM's{(t> = 0) eq.(40) we get in the a — )• limit 

ai oc — 61og(a/i), 03 oc — 31og(a/i). (122) 

However, the real task is to compute the constant terms here, which will give fj. 
The 3- and 4-loop ground state contributions are, analogously, 

Vr = + d,glXs + dsXl] 

Vr = - log(«/^)[ei53 + e2gt\s + esglXl + esA^], (123) 

where di and are pure numbers (we omitted for simplicity the contributions propor- 
tional to Xa)- Since there is (for a — )• 0) nothing to compensate for the dimensionalities 
of higher powers of coupling constants, there are no higher-loop contributions to Krac- 

6.6 The constants 03,01 and fj 

The constants 03 and ai can be found from the computation of the vacuum energy 2- 
loop diagrams (Fig. 1) in a theory with unbroken symmetry. This is equal to Kiac- The 
constant 7/ can be found from the 2-loop mass operator of the Higgs field containing the 
product of scalar and gauge coupling constants. In the Landau gauge the corresponding 
graphs are shown in Fig. 2. We denote this contribution to the mass operator as TiXg[p), 
where p is the momentum. It is easy to see that 

Sa.(0) = 3A3^. (124) 

From here we get a relation between and ?/, 

a3 = 3[log( — ) + 77 + c]. (125) 
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The constant ai is proportional to 03 with a simple symmetry coefficient, i.e. ai = 2a3. 
So, all three constants are related to one-another. We choose to compute the gauge 
invariant quantity Krac- This can be done in the unbroken phase, since this quantity is 
by definition ^-independent. The lattice gauge- Higgs vertices are shown in Fig. 3. 
The contributions of the two Feynman diagrams are 
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and ^ is the gauge-fixing parameter. The way to handle the trigonometric factors in 
the numerator is to separate from them terms that also appear in the denominator. 
Additional lattice corrections then also appear. Thus, for example. 
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and, putting for brevity a = 2, 

^ sin^(2pi + qi) = ^[2 sin^ Pi + 2 sin^(pi + qi) - sin^ gi - 4 sin^ pi sin^(pi + g^)] 

i i 

= 2{f + m') + 2(A;' + m') - g' - 4m' - 4 J] sin' sin'(pi + g^), (131) 



Z sin(2pi + qi) sin(2pj + g^) sin g^ sin g^ 

= E[sin^(^'' + - sin' Pi] [sin' (pj + g^) - sin'pj]. 



(132) 



^[sin'(pi + qi) - sin' Pi] 

= X][sin^ * + 2 sin Pi singi(l - 2sin'(^pi + ^gi))]. 



(133) 
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With these formulas and symmetries of the integrand one can firstly show that the 
term multiplying ^ in Isv + Issv vanishes: this checks the gauge independence of the 
result. In terms of the lattice integral I{m) one then obtains 

hv = Iglmmm) + ^{amfl{m)]. (134) 



In the MS scheme this vanishes since then 7(0) = 0. For the SSV diagram one has 

47(0)7(m) - P{m) - Am^HL{m,m,0) 
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where Hl is the lattice sunset function. The problem here is the last term; we find 
numerically, by expanding for small m, that it contains the term — 4m^-^/167r^, S = 1.94 
(see Appendix A. 2). Thus the total contribution to V^ac is 

Isv + Issv = 3gl[HL{m,m,0)m'' + ^7(0)7(m) + ^P{m) -6m'' + ^(am)'7(0)7(m)]. 

— . . . . ^^^^^ 

In the MS scheme this reproduces the glml term in eq.(40). Now the additional terms 

lead to the final result 

a3 = 3[log(^) + C + ^S^-^]. (137) 

This gives fj = 1.01. 

The most difficult task is the determination of the constant rj. Its computation with 
lattice perturbation theory needs quite complicated 3- and 4- gluon vertices. We do not 
attempt to make the corresponding computation. Instead, we determine this constant 
with sufficient accuracy in MC simulations. The idea of that computation is simple. 
One should pick up some quantity which can be perturbatively computed (at least to 
2-loop accuracy) and compare it with a result of Monte Carlo simulations. We choose 
the Higgs scalar condensate, defined below on the lattice, for these purposes. Its 2-loop 
continuum computation follows directly from the effective potential, while it can be 
determined by MC simulations with sufficiently high accuracy. 



6.7 The quadratic condensates on the lattice 

Any condensate defined in the MS scheme has its analogue on the lattice. The exact 
relationship between the lattice and MS condensates follows from the results of this 
and previous sections. We will give an explicit form for the quadratic scalar condensate 
only, the relation between higher condensates requires computation of the ultraviolet 
divergent vacuum graphs in the lattice regularization scheme on the 3- and 4-loop levels 
(i.e. a computation of the coefficients di and in eq.(123)). Using the exact relation 
between the lattice and continuum effective potential we get for the scalar condensate: 

S 3 6 
(</>V)iatt = 7^ + T^i^sClog — + fj + c)] + {<I>^<I>)r, (138) 
Z7ra lOTT"' aa 
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or, in terms of the radial mode of the lattice scalar field: 

Knowing the effective potential, one also can obtain a width of distribution of the order 
parameter on finite lattices. This point is discussed in Appendix C. 
For the Aq condensate we have, in complete analogy: 

{AIAIU, = 7^ + At^sClog — + fj + c)] + {AlADn, (140) 
47ra 47r"^ a/j, 

or, in terms of the lattice notations 

Of course, the lattice condensates are /is-independent, but divergent in the continuum 
limit a — )• oo. 

We have now defined the relation between (^V) in the continuum and on the lattice 
and computed it to 2 loops in the continuum. The comparison of the continuum 
computation with MC simulations will allow us to determine completely the constant 
physics curve. In fact, this comparison will provide information on the magnitude of 
higher order effects (3 loops, etc.) as well as on the magnitude of 2-loop finite scaling 
corrections (see below). 



7 Monte Carlo determination of the constant physics 
curves 

Now we are ready for a Monte Carlo determination of the parameter r] - the only one 
left in the relation of the lattice regularization scheme to the continuum one. 

Consider the theoretical prediction (139) for (-R|). This relation is exact in the 
continuum limit. In real MC simulations, one has a finite volume together with a finite 
lattice spacing a. So, there are finite size as well as finite a corrections to (139). Finite 
size corrections are easy to deal with, because in the broken phase with a non-zero 
mass gap m finite size corrections die exponentially with lattice size, ~ exp(— miV). 
Hence, for any fixed fia one can choose a volume large enough so that the value of (-R|) 
is volume independent. In other words, MC simulations can always be done in such a 
way that finite size corrections are not essential, and for fia < 40 the necessary volume 
is not very large. For the discussion of the probability distribution for the quantity 
{(f)^(f)) in the finite volume, see Appendix C. 

The finite a corrections are power-like and, therefore, more important. A way to 
improve the situation is to use the analytical lattice 1-loop expression for the average 
{<f)^<f)), generalizing eq.(77): 
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Aira 
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Then the unknown finite a effects appear on the 2-loop level only. With the use of 
eq.(142) the theoretical prediction for {(f)^(f)) is 

(</>V) = + (</>V)fi) + (</>V)(2) + ^(</>V), (143) 

where S{(f)^(f)) represents the value of the unknown higher order continuum contribu- 
tions. 

Let us analyse in more detail the behaviour of the function A = /^G/^j^Ai?^ = 
8/gl * A{(f)^(f)), where AR^ is the difference between the MC and the analytical results 
(143) (we take fSn corresponding to a "classical" regime, where spontaneous symmetry 
breaking appears on the tree level). There are three mass scales in our theory in the 
broken phase. Two of them (m^ and mi) are related to each other through ffiT jrhi = 
mw jfUH- The third one is the Debye screening mass mj). The general structure of A 
is: 

A = f[amT, ami)) + [gla)ip[ ; arriT, ami)), (144) 

TrrriT 

where /[arriT, ami)) represents the 2-loop finite a effects and ^ the 3- and higher-loop 
finite a effects. Here p = 5'|/(7rmr) is the expansion parameter of the 3d theory. The 
finite scaling behaviour of the function / is quite complicated due to the presence of 
two different mass scales, rriT and m/j, which have different dependence on the scalar 
mass ms. We notice, however, that if the vev of the Higgs field is large enough, then 
rriT — mj), and the function / essentially depends on one variable only. In the following 
we used the constraint v[T)/T > 3.6, which ensures that rriT — mj) with an accuracy 
of ~ 10%. In fact, this condition, being a bit arbitrary, has a weak influence on the 
parameter r] extracted from MC simulations. A further simplification comes about 
when we notice that for these values of vjT the loop expansion parameter is quite 
small, p ~ 0.1, and we neglect the 3- and higher-loop finite a effects. 
In general, there are three possibilities: 

(i) The lattice spacing is so small that finite a effects are not essential. Then the dif- 
ference between the lattice result and the 2-loop expression is related to the amplitude 
of the higher-order corrections to the continuum scalar condensate, 

A, = -^^(</>V), (145) 

93 

so that the function A2 depends on the variable gl/rriT = 2gT / (j) only. It must go to 
zero for gT / (j) — > 0. The value of rj can be found from this requirement. This type of 
behaviour is expected for a sufficiently small parameter am^, measuring the magnitude 
of finite a effects. 

(ii) Higher-order corrections are so small that their effects on the scalar condensate 
are negligible. Then the function A2 depends on the variable arriT only and must go 
to when amT/2 = (j)/[gTf3Q) — )• 0. The value of rj follows from this requirement. 
One expects to enter this regime when the expansion parameter p = gl/{T^mT) is small 
enough. 
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(iii) In the most general case one can neglect neither higher- order nor finite a correc- 
tions. Then the function A2 depends essentially on two variables, arriT and gl/rriT- 

We shall study the most general case. So, we need some expression for the finite a 
effects described by /{arriT) and an ansatz for higher-order corrections to the scalar 
condensate. The function f{x) can be expanded in a power series in a;, starting with a 
linear term: 

f{x) = Ax + Bx^ + .... (146) 

The linear term appears in the computation of the sunset diagram and the figure-of- 
eight diagram on the lattice (see Appendix A). Moreover, from a computation of the 
figure-of-eight graph we know that the (am)^ corrections come with a large coefficient, 
of the order of the coefficient in front of am. Therefore, in the analysis of the 2-loop 
finite a effects it may not be sufficient to consider the linear term only. So, we will 
keep two terms in the expansion of f{x). 

Simultaneously, we will use the following expression for the function 8[(j^(j)) account- 
ing for the 3- and 4- loop contributions: 

= -#7r^(/^ + 77^), (147) 

8A3 (47rj'*mr 47rmr 

where (3 and 7 are some unknown parameters. 

To summarize, the "theoretical" prediction for {R\) contains five numbers (?/, A,B, 
and (3 and 7) which we want to determine by comparing the prediction with lattice 
MC data. 

The parameters of the lattice action (83) are fixed as follows. We choose fia = 
12,20,24,32,40, thus fixing a (since g = 2/3) and mn = 80 GeV fixing A. The Aq 
couplings are fixed using eqs.(88,90) and eq.(92) with /2D = P2D = 0. The system 
then is simulated for various values of fSn {f^R is fixed by eq.(87)) for lattice sizes 
N = 12,16,24,32,40,48. An example of distributions in Rl for jSc = 12, mn = 80 
GeV and iV = 24 is shown in Fig. 4. They show a single peak deep in the broken 
phase, which develops into a 2-peak structure at the critical value fSffc- For each peak 
the value of (-R|) with error is computed. 

In Fig. 5 we present A as a function of (j)/[gTf3Q) = amT/2 for the best fit set of 
parameters for the Higgs mass niH = 80 GeV. The effect of variation of the different 
parameters of the fit is roughly as follows: the change of f3 and 7 changes the deviation 
of the points corresponding to different fSa from the single curve, parameters A and B 
fix the form of that curve, while the change in rj moves the curve up or down. The 
fitted parameters are: 

77 = 2.18(6), A = 5.09(25), B = -3.40(23), /3 = -162(22), (148) 

the parameter 7 characterizing the 4-loop contribution to the effective potential cannot 
be determined to any good accuracy by this fit; the numbers given above are stable 
with respect to variation of 7 within the limits —700 < 7 < 700. The quality of the fit 
is quite good, %^ = 38.9/40. In the computation of the continuum scalar condensate 
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we took fi3 = 1.6mr, the value for which 2-loop corrections to the effective potential 
are minimized in the vicinity of the minimum. 

The coefficient rj must be independent of the Higgs mass. To check this, we performed 
analogous simulations for a different Higgs mass, niH = 160 GeV, with the following 
result {i3g = 20, 32): 

7] = 2.44(29), A = 0.86(30), B = 0.67(19), /3 = -317(84), (149) 

with %^ = 11/9. Again, the parameter 7 remains undetermined. One can see that the 
values of all the parameters besides rj changed a lot, while rj remains constant within 
error bars. The value of /is, minimizing the 2-loop contribution for this value of the 
Higgs mass, is fis = 0.39mr. 

To check the stability of the result with respect to the constraint v/T > 3.6 we made 
a corresponding fit including the simulations with v/T ~ 1. As expected, the quality of 
the fit is worse than previously^^ (the function / for small v/T depends essentially on 
two variables rather than one), but the fit parameters are consistent with the numbers 
given above, i.e. 

77 = 2.17(4), A = 5.10(18), 5 = -3.43(18), /3 = -151(12), 7 = -98(83). 

(150) 

In this case the parameter 7 is also determined. 

It is instructive to compare the finite scaling contributions to A, which are deter- 
mined by this fit through the coefficients A and B with a typical 2-loop contribution 
~ {^j^ ~ 1 for niH = 80 GeV. These corrections are of the order of the 2-loop 
continuum contribution already at arriT ~ 0.5. 

The knowledge of the parameter f3 allows one to get an idea of the magnitude of 
the 3-loop corrections to the effective potential, at least at sufficiently large ^, so that 
rriT ~ m£). We parametrize the 3-loop corrections in this region as 

^3 = J^ydmT + Vr', (151) 

where V^^™^ is the piece of the 3-loop potential, which is singular in the limit m2 — > 0. 
It can be found from the equations of subsection 6.1, taking into account the condition 
that the 0{h^) contribution to the ground state is finite. We get: 

Now, to relate f3 and f3 we use the relation of the 3-loop contribution to the scalar 
condensate {(f)^(f)) through the 3-loop effective potential given in eq.(72). The difference 
between f3 and f3 appears from the existence of simply connected diagrams contributing 
to the condensate. We present here only the numerical result of this computation: 

/3~/3-113. (153) 



'The for this fit is 67/54 d.o.f, which gives a confidence level of 0.11. 
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This gives for the 3-loop contribution (3 = —49(22). This number is quite reasonable 

2 

and agrees with the expectation that the true loop expansion parameter is about 

[2]. _ _ _ . . 

We conclude this subsection by an estimate of the systematic uncertainty in the 
determination of r] due to the poor knowledge of the 2-loop corrections to the Debye 
screening mass^^. From eq.(1.53) the change of the effective scalar mass m\ due to the 
change of the Debye mass mj) and parameter r] is 

^-I = -A%-^^- (154) 
^ 167r2 16 ' 327rmz) ^ ' 

Using eq.(92) the variation of the Debye mass is 



Therefore, 



'5 



9l . , A 



o47r rriD 

So, for \pg + ^Pg\\ < 14 the systematic uncertainty in rj is the same as the statistical 
error. 



8 The theory with the Aq field integrated out 

In the main part of the paper we considered the 3d theory derived by dimensional 
reduction from the full 4d high-temperature one. In addition to the 3d Higgs field and 
gauge field it contains a triplet of scalar fields Aq. The mass of this field ~ gT is larger 
than the typical 3d scale ~ g'^T, so that this field can be integrated out perturbatively. 
This has been done in ref.[2]. The result of this integration is a 3d gauge- Higgs system 
with effective parameters related to that of the original 3d theory, containing Aq. We 
refer here to the relations (1.51,1.52,1.53) of [2]. 

In this section, we establish the constant physics curves for a corresponding 3d lattice 
gauge- Higgs system, i.e. we will relate it to the original 4d high-temperature one. 

The lattice action is given by (83), where all terms containing Aq are omitted. The 
continuum action is given by (1.51). The connection between the lattice couplings 
and the continuum parameters gs"^ and A3 is 

4 1 

f^G = (157) 
4 gi Pg 

^®There are no sizeable systematic uncertainties associated with the final volume effects, since in 
all cases simulations were done in such a volume that the final volume shift of {R\) was smaller than 
the statistical error. 



36 



In analogy with eq.(97) we write: 

mjj _1 1 
~ 2^° ^ {TaY 



3-^ + ^-ra/°" 



(3h (3h Stt 



fO m fO _ pO 

)^iog- + (ra) 32^^ 



(159) 



where 



lOTTi 

flm = + A3, (160) 



/2m = §53+9A353-12Al, 



Am = ^53^0 + 9A353'^ - 12A^77 + /^C 

Id 

Here the parameters fj and ?/ are the same as in the theory with the Aq field, but rjo 
is different from rj. The value of rjo can be defined in MC simulations in precisely the 
way we found rj in the previous section. We get: 

7/0 = 1.62(7), A = 2.62(14), 5 = -2.05(13), /3 = -67(20), (161) 

The fit (with = 24/25 d.o.f.) is shown on Fig. 6. We took /la = 2.37mr for 
these computations [2]. The parameter 7 remains undetermined, as in the previous 
cases. Now we can, with the use of eqs. (151,152), estimate the amplitude of the 
3-loop corrections to the effective potential in this theory. In complete analogy with 
the previous discussion we obtain 

/3 = /3- 52 (162) 

which gives f3 = —15(20). The error here can certainly be decreased by the increasing 
statistics. 

In fact, the relation between rj and rjo can be found analytically. To this end we 
computed the Aq contribution to the effective potential in lattice perturbation theory. 
The corresponding diagrams are shown in Fig. 7. The difference between the MS Aq 
2-loop contribution and the lattice one is found to be 

2 167r"' 16 a/is 5 4 8 

+ T#i6ml,[log( — ) + ^-S + Cl (163) 
167r"' a/is 4 

where the number p is related to a 2-loop integral defined in Appendix A and computed 
to be /) = —0.314. This relation, together with (1.53), gives 

81 51 ^,S2 15,, 3r 3 

—77 = —770 + 3( P-S) + —(log + — + ()■ (164) 

16 ' 16 ' M ^ ^ 8 ^ ^ 2mD 10 ^ ' 
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This establishes the required connection. With the parameter rjo determined above we 
can derive an independent estimate of rj: 

7] = 2.03(5). (165) 

This value is about 2 standard deviations below that defined in eq.(148). The difference 
between the two numbers may have an only statistical origin. Also, it may come from 
the systematic uncertainties associated with the Aq field (see the discussion at the end 
of the previous section). In particular, if pg + ^Pg\ ~ 15, the discrepancy disappears. 
In any case, the 2a uncertainty in the parameter rj is 0.12, which gives already 0.3% 
accuracy in the determination of the temperature through the lattice parameters. In 
absolute units, this means ST ~ 0.5 GeV in the vicinity of the phase transition for 
ms = 80 GeV. 

We conclude this section by noting that the 2-loop computations of the gauge in- 
variant condensates in this theory can be extracted from sections 5 and 6 of this paper 
by simply dropping the Aq contribution. 

9 Conclusion 

This paper is a quite technical (but absolutely necessary) step in the study of the 
electroweak phase transition with the use of the lattice MC methods in the framework 
of 3d effective theory. It provides a bridge between the lattice and continuum in terms 
of the constant physics curves. These curves are parametrized by 3 pure numbers 
{rj^fj and rj and are exact in the continuum limit. In addition, we found a number of 
relationships, which are exact in the continuum limit, between lattice and continuum 
gauge-invariant observables - condensates. The condensates were computed on the 
2-loop level. This gives a possibility to study the convergence of perturbation theory 
in the broken phase in the vicinity of the electroweak phase transition. 

The authors thank I. Montvay for many helpful discussions on dimensional reduction 
and on different aspects of lattice perturbation theory. K.F. is partially supported by a 
CEC program (CHRX - CT93 - 0319), K.R. is supported by United States Department 
of Energy grant DE-FG02-91ER40661. 

A Some 1- and 2-loop computations on the lattice 

In this appendix we derive the expressions for the 1-loop tadpole graph and for a 
number of 2-loop graphs on the lattice. 

A.l 1-loop graphs 

Abbreviating the lattice propagator by 

d{ni,n2,ns,m) = sin^(7rni/iV) + sin^(7rn2/-/V) + sin^(7rn3/iV) + (am/2)^, (166) 
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ni = l,...,iV — 1, we can define the lattice sum corresponding to a tadpole graph as 



1 1 



aI(N,m) = —-y— -. (167) 



In the limit iV — )• oo the sum can be converted into an integral, 

1 



1 /■"' 

aI(oo,m) = alim) = — - / d x- 

Att^ Jo i 



sin^ Xi + sin^ X2 + sin^ Xs + (am/2)^ 
= 1 ^cZae-H-)^[e-|«7o(la)]^ (168) 

4: Jo I 

where /q is the modified Bessel function. In the continuum limit a — )• 0, 

Km) = ^[S -am- ^{amf + 0.82^{amf + 0{{amY))] (169) 
47ra 



wi 



th 



S = -(18 + I2V2 - I0V3 - 7V6)K((2 - V^YiVs - V2Y) = 3.1759114, (170) 

TT 

K being the complete elliptic integral of the first kind, and ^ = 0.15281. The values of 
^ and of the next coefficient are the results of a numerical computation. 
A related function is the 1-loop contribution to the vacuum energy, 

1 1 

a^J{N,m) = ^^og[d{ni,n2,ns,m)]. (171) 

ni=0 

In the continuum limit J(oo, m) = J[m) satisfies 2dJ j dm^ = I{m'^) and has the small-a 
expansion 

J(m) = const + f-Sm^ - -am^ - -^aV^ + 0{a^m^)\ . (172) 
47ra \2 3 4 / 

A. 2 2-loop graphs 

The scalar sunset integral H is given in eq.(1.22). Its lattice analogue is the double 
sum 

-j^ N-l N-1 

H{N,ml,m2,m3) = 7^ E E W^i'^2,n3,ml)]-i (173) 

nj=0 mj=0 

[(Z(mi,m2,m3,m2)]~-^[(Z(ni + mi, ^2 + m2, ^3 + m3,m3)]~\ 

where < + < iV — 1. The iV — )• 00 lattice limit [N — )• 00, a — )• 0, Na = 
constant) is also given by 

Hi,{rn,m,m) = if (00, m,m,m) 
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^ r-^ 1 1 1 

— — / (fx(fy = = = 

647r^ Jo Yl sin^ Xi + ^a^m^ Yl sin^ Vi + ^a^m^ Yl sin^(a;i + yi) + ^a^m^ 

1 1'°° 1 

- / d^a exp[-3(l + -a^m^)(ai + a2 + "s)] • 
o ^0 Z 



8 

,2,r (^a; , 1 3 



(174) 



/* dx I 

/ — e"^"°'^7o(va? + ai + Saiascosa;) 
Jo 27r 

~ flog — + ^ + C + Ka"^) + C'(M')) 

^ — + log ^ + ^ + C + + 0{{amf)\ . 

Vxi-K^ \ afis dm 2 J 

The constants ( ~ 0.09 and u ~ —0.6 have been estimated by computing the integral 
(174) numerically and also by computing the sums for values of N up to 24. Since 
the dependence on am/2 is only logarithmic, one has to go down to am/2 ~ 0.01 
to separate the logarithmic and constant terms, and the computation is numerically 
demanding. 

We finish this appendix by presenting two non-trivial 2-loop integrals needed for the 
computation of the parameter rj and for relating rj to rjo (the parameter for the theory 
with the Aq field integrated out). Both come from the SSV type of graphs. 

The number 6 is related to the derivative of the SSV diagram with respect to the 
scalar mass and is given by 

r ^ f' js f' js E»[sln^ Xj sin^(a; + y)^] 

" = 7r~r / ax a y—— — r-^ — tttzz — t-^t ttt^^ — r-^ — r — 1-94. (175) 

27r4y_| y_| '{EiSm'xiy{YiSm\x + y)i){YiSm'yi) ^ ' 

The number p is related to the derivative of the SSV diagram with respect to the 
vector mass and is given by 



^ J_ /-f ^3_^ /-f ^3 f E»[sin^a^»sin^(a; + y) 
^ " 47r4 y_| y_| ^ \ {Yi sin' Xi){Yi sin'(a; + 



E- 4 
•sm Xi 



y)i) (Ei sin' a;i)2 j (Eisin'yi)2 
0.314. (176) 



B The gauge invariant effective potential 

A gauge-invariant effective potential V{(t) has been introduced in [13] by including a 
current J coupling to cf^cf): 

exp[-ir(J)] = j VA^Vcf) exp[-5 + j d^xJcf^^] (177) 

and by performing the usual transformation from J to a. Thus simply (the three- 
volume is not written explicitly) 

W{J) = V{v{T\ml + J\ (178) 
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i.e. the generating functional W{J) is the same as the effective potential at its mini- 
mum, but computed for mg — )• ml + J. The relation between a and J is given by 

WV) - 5I(^M±i) . ^ . (179) 

and the gauge-invariant effective potential is given by 

V{a) = V{v{T),ml + J{a)) - J{a)a. (180) 

To evaluate the Legendre transformation (180) we write J = J(o) + {J{2) is not 

needed). To leading order in it follows from eqs.(71) and (179) that 

J(o) = -ml - 2A30- (181) 

so that from eq.(71) 

V(o)(ct) = V{v{T),ml + J(o)) = mla + \^a\ (182) 

To order E 



where now 



3^1 _ ; 

mi + — (2mr + mi) 
0A3 



(183) 



'1 „ „ 1 



mi = ■^4A30-, mr = y -^-lo", mL = ^m% + -gla. (184) 

These equations show that the calculation is valid for <t > or for ml + J < 0, i.e. in 
the broken phase. With these masses, the 2-loop result for the gauge invariant effective 
potential is 

V{a>0) = mgo- + Aso-^ - -|-(6mr + 3m| + mi) (185) 

127r 

-— — mi[5mi + — -^(2mr + m^)] + ti^V2{fhT,fhL,fhi,m2 = 0), 
D47r"' zAs 

where V2 is given in eq.(34). 

If <T < or ml + J > one has to use the eqs. (73), (74) and (40), which are 
appropriate for the symmetric saddle point. Now the leading term in V{v) in eq.(40) 
is 0{h) and one obtains 

h 



<T{J) = - — ^^ml + J. (186) 

ZTT 

The same procedure as before leads to 

V{a < 0) = (ml - ^^glmn) a + ^ (^3^ log ^ + \gl + 2A3) - (187) 

neglecting a cr-independent part depending on m/j, log(/i3/m/)), and the coupling con- 
stants. Note the dependence on h. 

Although explicit forms of the gauge-invariant potential have been given here one 
has not gained anything of practical use. One is anyway only interested in the physical 
ground-state value of {(f)^(f)) and this is already given by eqs. (76-78). 
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C Distribution of ^ finite system 

In the course of lattice simulations = J (Px(f)^[x)(f)[x) will be studied in finite 
volume systems and it is of great use to know its distribution dN j dB? in the continuum, 
but in a finite 3-volume V3. By definition we have, in terms of the functional integral 
in eq.(38). 



dN 



= r ^ f'^i' exp(-5[^] + itR^ - it^ [ d^x(t>^) 
J -00 2ir J Vs J 

/oo M 
— exp{itR^ - VsV[v{T),ml + it/Vs]}. (188) 
-00 2Tr 

The integral over t is computed with the saddle point method: the saddle point ts is 
determined by 

dV[v{nml + itjVs] dV[v{nml] .t, d'V[v{nml] 

dml ~ dml ^Vs {dmlY ^ - ' ^ ^ 



so that the saddle point is at 



it, R^ - V'[ml] 



Vs V"[ml] 



(190) 



where the derivatives always are with respect to mg, keeping (j) = v[T). Near the 
average 

{R') = V'iml), (191) 

ts/Vs is small relative to ml and it is consistent to expand as in eq.(189). For the 
distribution one obtains 

= =exp| ^ IR' - V'imDA. (192) 



This is a Gaussian with the average as in eq.(191) and the width —V"{ml)/V3. Note 
that V"{ml) < 0: for a simple tree potential near a broken minimum [ml < 0) 
the derivatives are V{ml) = -m^/(4A3), V'{ml) = -ml/{2X3), V"{ml) = -1/(2X3). 
Eq.(192) shows how the fluctuations freeze to zero in the infinite volume limit. 
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Figure 1: The vacuum 2-loop diagrams used in the computation of r\. 
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gure 2: The Higgs self-energy 2-loop diagrams used in the computation of 
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6ij cos la{pi +p^) 



Figure 3: The scalar-gauge vertices in lattice perturbation theory. 
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Figure 4: The evolution of the distribution of (-R|) with fSn- 
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Figure 5: The quantity as a function oi x = amT/2. 
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Figure 6: The same as Fig. 5, for a theory without Aq field. 
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Figure 7: The Aq contribution to the effective potential. 
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